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Abstract 

By making use of the path-integral duahty transformation, string representation of the Abehan- 
projected S'f/(3)-QCD with the 9-term is derived. Besides the short-range (self-) interactions of 
quarks (which due to the G-term acquire a nonvanishing magnetic charge, i.e. become dyons) and 
electric Abrikosov- Nielsen- Olesen strings, the resulting effective action contains also a long-range 
topological interaction of dyons with strings. This interaction, which has the form of the 4D Gauss 
linking number of the trajectory of a dyon with the world-sheet of a closed string, is shown to 
become nontrivial at B 7^ Svr x integer. At these values of 0, closed electric Abrikosov-Nielsen- 
Olesen strings in the model under study can be viewed as solenoids scattering dyons, which is the 
4D analogue of the Aharonov-Bohm effect. 

PACS: 12.38. Aw; 12.38.Lg; 11.15.-q 

Keywords: Quantum chromodynamics; Effective action; Confinement; Aharonov-Bohm effect 



*Permanent address: Institute of Theoretical and Experimental Physics, B. Cheremushkinskaya 25, RU-117 218 
Moscow, Russia. 

^E-mail address: antonov@mailbox.difi.unipi.it 



1 



1 Introduction 



Recently, by making use of the method of Abehan projections a remarkable progress in the 
description of confinement in QCD has been achieved (for a recent review see [|]). In this way, 
most of the results have been obtained for the simplest case of the S'[/(2)-group. That is because in 
this case the respective Cartan subgroup is just U{1). Owing to this fact, under the assumptions 
on the irrelevance of the off-diagonal degrees of freedom in the IR limit (the so-called Abelian 
dominance hypothesis) and condensation of monopoles arising during the Abelian projection, 
the resulting effective theory takes the form of the dual Abelian Higgs model. In the latter model, 
confinement of external electrically charged particles is provided by the formation of electric 
Abrikosov-Nielsel-Olesen (ANO) strings 0] between these particles. This is just the essence of 
the 't Hooft-Mandelstam scenario of confinement |^. 

Besides this simplest S't/(2)-case, an effective [f/(l)]^ gauge invariant Abelian-projected the- 
ory of the realistic S't/(3)-QCD, based on the Abelian dominance hypothesis, has also been pro- 
posed [^]. Quite recently, confining properties of this model were investigated by a derivation of 
its string representation ^ (for a review see Moreover, in Ref. the collective effects in 
the grand canonical ensemble of strings in this model were investigated. In particular, a crucial 
difference of this ensemble from the one of the SU{2) Abelian-projected theory |jlO| (following 
from the appearance of two types of (self-)interacting strings in the SU (3)-case) has been found. 

The present Letter is devoted to the investigation of Abelian-projected SU (3)-QCD extended 
by the introduction of the B-term. The relevance of this term to the instanton physics is a well 
known issue, and we will not discuss it further, referring the reader to classical monographs (see 
e.g. |n[]). Within our investigations, it will be demonstrated that the introduction of the G-term 
leads to the appearance of a long-range topological interaction of quarks (which due to this term 
become dyons by acquiring a nonvanishing magnetic charge |T2|) with the world-sheets of closed 
electric ANO strings. This interaction has the form of the Gauss linking number of these objects 
and can be interpreted as the 4D analogue [jl3| of the Aharonov-Bohm effect where strings 



play the role of solenoids, which scatter dyons. In particular, it will be demonstrated that at 
certain discrete critical values of 9 this effect disappears. This can be interpreted in such a way 
that at these values of 0, the relation between the electric flux in the ANO string (viewed as a 
solenoid) and the magnetic charge of the dyon is so, that the scattering of dyons over strings is 
absent. Note that the long-range interaction of dyons with magnetic ANO strings in the Abelian 
Higgs model (which corresponds to the theory dual to the Abelian-projected SU{2)-QCD with 
the 0-term) has been studied in Ref. |T^. Similar interaction of an external electrically charged 
particle with magnetic ANO strings in the same model but without B-term has been found in 
Ref. |]l6l by the evaluation of the Wilson loop describing this particle. 

It will be also shown that, as it can be expected from the very beginning, besides this long- 
range interaction of strings and dyons, the resulting effective action contains the short-range 
(self-)interactions of these objects. The short-rangeness of all these interactions is due to the 
(self-) exchanges of strings and dyons by massive dual gauge bosons. Among these short-range 
interactions, the selfinteraction of the open world-sheets of strings, which end up at dyons, is the 
most important one, since it is this interaction, which provides the linear growth of the confining 
dyon-antidyon potential. 

The organization of the Letter is as follows. In the next Section, we shall discuss an effective 
theory of the Abelian-projected S'?7(3)-QCD with the O-term. In Section 3, a derivation of the 
announced string representation for the partition function of this model will be presented, and the 
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obtained results will be briefly discussed. Finally, in three Appendices, some technical details of 
calculations, which are necessary for a derivation of the formulae from the main text, are outlined. 



2 The Model 

In the present Section, we shall discuss the model, whose string representation will be derived 
below. To start with, let us consider the partition function of an effective [[/(l)]^ gauge invariant 
Abelian-projected theory of the pure S'f/(3)-gluodynamics {i.e. a theory without quarks and the 
e-term) ||. It reads Q 

r / \ 



a=l 

2X2 



\a=\ 

X exp |- 1 ^Fj, + [K^M - ^Sm^a^t.) ^af + A (|<I>„P - 1] 

Here, $a = |$a| e*^" are the dual Higgs fields, which describe the condensates of monopoles emerg- 
ing after the Abelian projection, and gm is the magnetic coupling constant related to the electric 
one via the topological quantization condition = 4:iTk. In what follows, we shall for simplicity 
restrict ourselves to the monopoles possessing the minimal charge only, i.e. set k = 1. Next, 
is the "magnetic" potential dual to the "electric" potential = (^A^, , where A^^ are 
just the diagonal components of the gluonic field. On the R.H.S. of Eq. (p, there have also been 
introduced the so-called root vectors e^'s, which have the form 





ei = (1,0) , 62 = 

These vectors play the role of the structural constants in the first of the following commutation 
relations 

[H, E±a] = ±eaE±a, [E±a, E±b] = ^-j^SabcE^c, [Ea, -^-b] = ^ab&a}^- 

In these relations, H = {Hi,H2) = (T^,T^) are the diagonal S'f/(3)-generators, which generate 
the Cartan subalgebra, where from now on T* = y, i = 1, . . . , 8, are just the S't/(3)-generators 
with A"s denoting the Gell-Mann matrices. We have also introduced the so-called step operators 
E±aS (else called raising operators for positive a's and lowering operators otherwise) by redefining 
the rest (non-diagonal) 5'f/(3)-generators as follows 

(Clearly, these operators are non-Hermitean in the sense that {EaY = E^a-) It is also worth 
noting that owing to the fact that the original SU{?>) group is special, the phases 6'a's of the dual 
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Higgs fields are related to each other by the constraint I] 6'^ = 0. This constraint was imposed 

a=l 

by introducing the corresponding 5-function into the R.H.S. of Eq. (|l|). 

^ Throughout the present Letter, all the investigations will be performed in the Euclidean space-time. 
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In what follows, we shall be interested in the study of the model in the London limit, i.e. the 
limit of infinitely large coupling constant A of the dual Higgs fields (Clearly, this limit corresponds 
to the case of infinitely large masses of the dual Higgs bosons.). It is just this limit, where the 
model under study allows for an exact reformulation in terms of the integral over closed electric 
ANO strings 0. In the London limit, the radial parts of the Higgs fields can be integrated out, 
and the partition function (|I|) takes the form 



Z 




X exp< / d^x 



^ 1 1 



(2) 



Since in the model (|l]) there exist string-like singularities of the ANO type, in Eq. (|^) we have 
decomposed the total phases of the dual Higgs fields into multivalued and singlevalued (else oftenly 
called singular and regular, respectively) parts, 9a = 6'^™^' + 0^^^', and imposed the constraint of 
vanishing of the sum of regular parts by introducing the integration over the Lagrange multiplier 
k{x). Analogously to the dual Abelian Higgs model, in the model (0), 6'^™s 's describe a certain 
electric string configuration and are related to the closed world-sheets S^'s of strings of three types 
via the equation 



e 



HvXp'-'X'-'p 



X) 



2vrS;:,(a;) 



= 27: da 



(x('^)(0)5(x-x('^)(0 



(3) 



This equation is nothing else, but the covariant formulation of the 4D analogue of the Stokes 
theorem for the gradient of the field 9a, written in the local form. In Eq. (^, x^"''^{^) = x|f^(^) is a 
vector, which parametrizes the world-sheet with ^ = (^^, ^^) standing for the 2D coordinate. On 
the other hand, the regular parts of the phases, 9^a^''s, describe a singlevalued fluctuation around 
the above mentioned given string configuration. Note that owing to the one-to-one correspondence 
between ^^™s-'s and S^'s, established by Eq. (|^), the integration over 6'^™^ 's is implied in the sense 
of a certain prescription of the summation over string world-sheets. One of the possible concrete 
forms of such a prescription, corresponding to the summation over the grand canonical ensemble 
of virtual pairs of strings with opposite winding numbers, has been considered in Refs. 0, |10[- It 
is also worth noting that due to Eq. (^) the integration measure over the full phases ^^'s can be 
shown to become factorized into the product of measures over the fields ^^'°^''s and 6'Jj°S''s. 

As it was announced in the Introduction, the model whose string representation will be of our 
interest is not the London limit of pure Abelian-projected 5'f/(3)-gluodynamics (0), but rather 
the full SU {3)-QCD extended by the 0-term. For a certain quark colour c = R, B, G (red, blue, 
green, respectively), the partition function of such a theory reads B 



Z, = J VB^V9:-^-V9l'^-VkS (E ^f'-) expj J d'x ~ (f^, + F 



p,u 



E i9,9a - 9meaB,f + ikj^ 



a=l 



a=l 



16vr2 V*'^" 



p(c)^ (-p 



F(c) 



(4) 



^Note that the size of the core of the string (vortex in 3D) is equal to the inverse mass of the dual Higgs fields, 
which means that the London limit corresponds to infinitely thin strings. 
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Here, O^u = ^£fMuXpO\p, and F^^^ is the field strength tensor of a test quark of the colour c, which 



moves along a certain contour C. This field strength tensor obeys the equation d^F^^} = Q 

where = f (ia;^(r)5(x — x{t)), and the vectors of colour charges read 

c 



Q 



Clearly, F^fj = Q^^^T^u, where either 






47r2 



_d_ 

dx. 



A- 



jujy) _ _d_ 
X — nY dxy 



x-yf 



(5) 



or J^,, 



-Unu. Here, S, 



— —^piv ^iiv\^) = J d(7^iy{x{^))5{x — x{^)) is the vorticity tensor current associated 

E 

with the world-sheet S of a certain open electric ANO string, bounded by the contour C. 



3 Aharonov-Bohm Effect and Confinement 

In order to proceed with the string representation of the model @, it is useful to employ the 
so-called path-integral duality transformation (see e.g. Ref. |]T^, where this transformation is 
presented for the simplest case of the (dual) Abelian Higgs model), whose details are outlined in 
Appendix A. The result of this procedure reads Q 



exp 



127r2 



d^xj^L-^L{j:,c) 



vx^'^\o^J:^ 



X 



\a=l 



(6) 



Here, m. 



^fj-^ix + ^xhlj^u + ^i^^X/i is the field strength tensor of the Kalb-Ramond field of the 
a-th type, h"^^, J^^^ is defined by Eq. d), and L(S,C) = ^e^^xp I d^xd^yj:f,^{x)jp{y)-£^j^^ is 
the (formal expression for the) 4D Gauss linking number of the surface S with its boundary C, 
which will be shown to become cancelled from the final expression for the partition function. In 



Eq. 

numbers s^^^^s: s 
obey the relation Q 



), we have also denoted S^^ 

(R) _ AB) _ (G) _ 
3 — ''2 — ~ 
(c) = fip Jc) 



ya 



0, S\ 



(R) 



47r2 

■^2 



s^f'Tn.v and introduced the following 



'^2 



— S 



— S 



(G) 



1, which 



As the next step, one needs to carry out the integration over the Kalb-Ramond fields /i^^^'s. 



Referring the reader for details of this procedure to Appendix B, we shall present here only the 
outcome of the respective calculation, which has the form 



m;:,exp^ I d^'x 



24:7] 



exp 1-271^ J d^xdSvl^\x - y) 



r/2s;^(x)E;:,(y) + -— j,(x)j,(y) 



(7) 



•^Froin now on, we omit everywhere the normalization constant, implying for every colour c the normalization 
condition Zc[C — 0] — 1. 
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Here, 



X) 



m Ki{m\x\) 



47r^ 



is the propagator of a dual vector boson of the mass m = gmvV^ with 



Ki standing for the modified Bessel function. 

Further, it is necessary to carry out the integral 



emerging among other terms in the expression —2t:'^ rf j d^xd'^y'E^j^{x)'D!^\x — which 
stands in the argument of the exponent on the R.H.S. of Eq. (|^. This calculation is outlined in 
Appendix C, and the result reads 



1 fQg. 



6 V 27r2 



d xd yjf,{x) 



{x - yf 



An'vl^\x-y) 



(9) 



On the other hand, after substituting Eq. into the term — -^^fp^ / d'^xJ^'^^, which stands 
in the first exponent on the R.H.S. of Eq. (|^), we get for this term the expression 
— '•^4™i / d^xd'^yj^{x) (^^}yy j^{y), which precisely cancels the corresponding Coulomb interaction 
of currents in Eq. (^. Clearly, this cancellation could be anticipated from the very beginning just 
from physical principles, since in the resulting theory there does not remain any massless particles, 
by which the quarks might exchange. 

Substituting now Eqs. (J^) and into Eq. (||), we arrive at the following intermediate result 
for the partition function 



Zc = exp 



4ie 



67r2 



d^xd^yj^{x)V^^\x - y)j^,{y) 



X 



X 



J Vx^^\OS {T^^^Uyw ~2inr^yJ d'xd'y±%{x)Vl^\x - y)±;M + 

+ zQgiv'sl:^ J d'xd'ytl,{x)V^^\x - y)^,M , (10) 

where it has been denoted S^j^ = S^^^ — s^f'T^^j,,,. Clearly, the last term in the argument of the 
second exponent on the R.H.S. of Eq. (|l^) requires further simplifications. By virtue of Eq. @ 
and a partial integration, this term can be rewritten as 



2 Jl „(c) 



(47r 



2\2 



sf^ep.vxp I d^xd^yd^zi:'' (x) 



jxjz) d Ki{m\x-y\ 
{y — zy dxp \x — y\ 



The integral over y has the same form as the integral /(A) from Appendix C with \ = x — z and 
47r^X>^)(x — z) . Owing to this result and the definition of , the term under 



reads ^ 



^ (x—z)^ 

study takes the form 



2si^)L(S", C) - 4L(S, C) - s^h,,xp j d'xd%Vl^\x - y)j,{x) 



_d_ 

dy. 



One can now see that, as it was expected, the singular term — ^L(Z1, C) from the first exponent 
on the R.H.S. of Eq. ([T0|) becomes cancelled by the corresponding term from Eq. (|ll|). Finally, 
we arrive at the following expression for the partition function 
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exp 




\a=l 



X exp 



(c) 



-2{7rvf J d'xd'y%{x)V^}{x - y)tlM + 



d'xd'yV^^\x-y)3, 



[X 



_d_ 



2iQ 



(12) 



which is the main result of the present Letter. Notice that for every colour c, it is straightforward 
to integrate one of the world-sheets S'^'s out of Eq. (|1^) by resolving the constraint imposed by the 
corresponding ^-function. This procedure leads to an obvious expression for the partition function 
in terms of two independent string world-sheets. 

The first exponent on the R.H.S. of Eq. (|12]) represents a short-ranged interaction of quarks, 
which due to the B-term have now acquired also a nonvanishing magnetic charge, i.e. became 
dyons. Clearly, the first term in the second exponent on the R.H.S. of Eq. ([T^ is the short-ranged 
(self-)interaction of four different world-sheets: three closed ones S"'s and an open one S. This 
term is responsible for the linearly rising potential, which confines dyons. Upon the expansion in 
powers of the derivatives w.r.t. ^°'s, this term yields the coupling constants of the local string 
effective action (see the first paper from Ref. [0 for details). The second term is just the short- 
range interaction of dyons and strings. The most nontrivial term in the obtained expression is the 
last, third one, which describes a long-range interaction of dyons with closed world-sheets. Such 
an interaction represents the 4D-analogue of the Aharonov-Bohm effect [jTSf, which means that 
at 7^ Svr X integer (c/. the explicit form of the numbers s^'^^'s), the dyons become scattered by 
the closed electric ANO strings. Contrary to that, these critical values of correspond to such 
a relation between the magnetic charge of a dyon and the electric flux inside the string when the 
scattering is absent. 
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5 Appendices 



A. Path-Integral Duality Transformation 

In the present Appendix, we shall outline some details of a derivation of Eq. (P) from the main 
text. To proceed with, we shall firstly rewrite a 6'^'^^'-dependent part of the statistical weight, 
entering the partition function (El), as follows 



j pc^ mexpjy d' 



X 



a=l 
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a=l 



and carry out the integration over ^^''^-'g^ ^^j-^jg way, one needs to solve the equation d^C^ 
which should hold for an arbitrary index a. The solution to this equation reads 



d 



477^ dx. 



d'y- 



k{y) 
X - yf 



where /i^^ stands for the Kalb-Ramond field of the a-th type. 

Secondly, replacing the integrals over ^^™^''s with the integrals over x^"''^ (O's by virtue of Eq. 



and discarding for simplicity the Jacobians |jT6[ emerging during such changes of the integration 
variables 0, we arrive at the following representation for the 6'a-dependent part of the partition 
function iW) 



3 



\a=l 



a=l 



a=l 



Vkexpl^-^ j (fxd'^y 



3 k{x)k{y) . I d k{y) 



Arf{x-yy \dx^{x - yy J 



x] 



X 



X J (05 E J Vh;, exp J d'x 



\a=l 



247] 



'■ (Al) 



Here, H^^^ = d^h^y^ + dxh^^^^^ + dyh'l^ stands for the field strength tensor of the Kalb-Ramond field 

3 

/i"^. Clearly, due to the explicit form of the root vectors, the sum e^B^ vanishes, and the inte- 

a=l 

gration over the Lagrange multiplier k thus yields an inessential (field- independent) determinant 

3 

factor. Notice also that due to Eq. (|^), the constraint Q^'^s- = q j-^as gone over into a constraint 

a=l 

3 

for the world-sheets of closed electric ANO strings of three types J2 '^uiy = 0- This means that 

a=l 

actually only the world-sheets of two types are independent of each other, whereas the third one 
is unambiguously fixed by the demand that the above constraint holds. 

Let us now turn ourselves to the pure gauge field sector of the partition function (^. In 
this way, firstly the term — | / (i^xF^j,F|fj can be rewritten as Q*^"^^ / d^xB^duTj^y. Secondly, the 

9-term is equal to ^^^s^^xp^ d'^xB^^d^F'^j^. (It is worth noting that this expression can further 

be written as — ^^^'^ q^^^ / rf^xB^j^ with q*^^^ = being just a (^-independent vector. The last 
equation means that due to the B-term the quarls acquire a nonvanishing magnetic charge, i.e. 
become dyons. The G-term then describes an interaction of a dyon with the dual gauge field WM.) 
After that, the gauge field sector takes the form 



Alternatively, these Jacobians can be referred to the integration measures 
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JvB^exp'.J d'x -j(f„„ 



167r2 



"DB „ exp {— d^x 



Q 



(A2) 

It is further convenient to pass from the integration over the B^-fields to the integration over 
the fields i?^ = GaB^. Then, performing the rescahng i^^Cfi^^) = y|iJ|^(oi<i) order to restore 

the standard factor 1/4 in front of (^F^u) ) taking into account that Q'^'^^ = ^Gas'f^ (see the 
notations after Eq. ©), we get for Eq. (A. 2) the following expression 



+ 



P5;exp|-/rf^.[i(F;,)%l(F(r)) 



After that, the integration over the i?"-fields can be carried out by making use of the first-order 



formalism, i.e. by representing the factor exp 
antisymmetric tensor field as follows 



as an integral over an auxiliary 



exp 



^ pv 



The S^-integration then yields the following expression for the field 



pu 



ym"'pu ^ 2 ''a '^pu 



(c)f 

37r ^\ 



(A.3) 



with A^'s standing for the new fields, resulting from the resolution of the corresponding constraint. 

Bringing now together Eqs. (A.l) and (A. 3), we arrive at the following expression for the full 
partition function 

Z, = j Vx^^\i)5 (jz S^.) Vhl^VX^e^Y^l j d'x\-^^ (h;^,)' - 



Ka=l 



Qgm 

37r 



'/7 ; 



dpK 



- 2 

+ mKJll^ - - fF|;j) 



Clearly, due to the closeness of the world-sheets S°'s, the obtained action standing in the exponent 
on the R.H.S. of Eq. (A. 4) is hypergauge- invariant, i.e. it is invariant under the transformations 
/il, — * h^^i, + 5^7" — d^j^, ^ A^ — ^\f^l^- Owing to that, the fields A^'s can be completely 



eliminated by setting 7^ 



|A^, which leads to the following expression 
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Next, a straightforward algebra yields 



(A5) 



3n 



where we have taken into account that for every c, (s^-^)) = 2. Clearly, the sin- 
gular term — | (F^fj^ becomes cancelled by the term fg^^^- Finally, substituting for 
J^^i, Eq. (H), we can rewrite the term ^ J d^xT,^^^^^ as — ^I/(S,C), where L(T,,C) = 
-^SfiuXp J d'^xd'^y'£^^{x)jp{y)-£^ ^^}y^2 is the formal expression for the 4D Gauss linking number of 
the surface S with its boundary C. As it will be demonstrated below, this singular term eventually 
cancels out. 

Accounting for the above considerations in Eq. (A. 5), we arrive at Eq. of the main text. 



B. Integration over the Kalb-Ramond Fields 

In this Appendix, we shall present some details of a derivation of Eq. (|^ from the main text. 
Namely, we shall carry out the following functional integral over the Kalb-Ramond fields 



Vh^exp 



d x\ 



Tja 



+ 



3^, 



{B.l) 



Clearly, such a Gaussian integral can be calculated upon the substitution of its saddle-point value 
back into the action. The saddle-point equation in the momentum representation reads 



This equation can most easily be solved by rewriting it in the following way 

where m = gmvV^ is just the mass of the dual gauge bosons, equal to the mass of the Kalb- 
Ramond fields. In Eq. (B.2), we have also introduced the following projection operators 

P^iu,\p = 2 i'Ppx'Pup — 'PfjtpVux) and lp,v,xp = 2 {^p.x^yp ~ ^p.p^ux) 
with V^u = Spu — These projection operators obey the following relations 
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^lj,u,Xp ^ufj.,Xp l/^i/,pA 

(the same relations hold for Pfj,u,\p), and 



^pu,af) 



P. 



p.v^\p 



ii-P) 



Xp,af3 



0. 



(5.3) 



(5.4) 



By virtue of the properties (B.3) and (B.4), the solution to the saddle-point equation (B.2) reads 



1 + 



1-P 



- Xu,al3 



which, once being substituted back into the functional integral (B.l), yields for it the following 
expression 



exp< —2T:'^rf' 



1 + 



1 - P 



af3\ 



(5.5) 



- pv,af3 



Rewriting now the term from Eq. (B.5), which is proportional to the projection operator 1, in the 
coordinate representation we obtain the expression 2n'^rf ^ d'^xd'^yT,'^^j^{x)Vl^ {x — ?/)S^^,(?/) from 
the action standing in the exponent on the R.H.S. of Eq. (|^). 

Let us now consider the term proportional to the projection operator — on the R.H.S. 
of Eq. (B.5). In this way, by making use of the following equation 

P^{i- P)fiu,aP = \{^upPpPa + Sf,aPuPl3 - ^uaPpPp " ^^^PuPa), 



we see that the term under study. 



2 2 / d^P 



P 



(27r)^ p^ + m? m? 
after carrying out the integration over p, reads 



(i - p) / rfSdV^(^-^)s^,(x)s^^(y), 



27irj\ 



m 



) ^ / d'xd'yf:;, (x) (y) Q^T^t^ i^-y). 



Performing twice the partial integration and making use of the equation = s^^'j^, we can 

further rewrite this expression as — |^ / d'^xd'^yj ^{x)V':^ {x — y)j^{y), which is just the last term 
in the exponent on the R.H.S. of Eq. (0). 

C. Calculation of the Integral (|^) 

Below in this Appendix, we shall present some details of a calculation of the integral (^. Owing 
to Eq. (|^) and the fact that for every c, (^s'f''^ = 2, this integral can be rewritten as follows 



d u- 



Stt^ \^ 47r2 J J ^ k — 1/1 dx^ J {x — zY dy^ J {y — u)"^ 

Next, by making use of the following sequence of partial integrations 



(C.l) 
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(9 1 d Ki{m\x — y\) d Ki{m\x — y\) d 



dy^ {y - uY dyf, \x - y\ dx^ \x - y\ dx^ {x - zf ' 

one can carry out the integration over z. After that, the integral (C.l) takes the form 



27r2 \ 47r2 



d'^xd^ujf,{x)I{x - u)j^{u), (C2) 



where I(x — u) 

= / ^^y \l-^{y^$ is just the integral, which is left to be calculated. 
In order to proceed with this calculation, let us pass to the integration over the variable 
z = y — u and introduce the vector X = x — u, after which we obtain /(A) = / g;4^ -f^i(H^-A|) 
Such an integral can most easily be carried out by substituting for the two ratios, which form the 
integrand, the R.H.S.'s of the following equalities 

Ki(m\z-\\) 47r2 r d^p ^p^^-^) 1 o f d^q ^''^^ 

- I ^ and — = 47r^ ' 



Then we have 

/(A) 



A| m J (27r)4p2 + ^2 ^2 J ^27t)^ 

IQtt^ f d^p e-'P^ 



m J (27r)^ (p2 -)- rn?)p'^ 



167r^ 



hoo +00 



00 +00 -/3m^- ,/, 

g ' 4(c«+/3) 



m J (27r) 



(C.3) 



It is further convenient to introduce new integration variables a G [0, +00) and t G [0, 1] according 
to the formulae a = at and P = a(l — t). Straightforward integration over t then reduces the 
integral (C.3) to the following expression 

r, +00 



Such an integral can be carried out by virtue of the formula 

+00 ' 



x'^-^e-^-^'^dx 



and the result has the form 



/(A) 



.7, 



47r2 
m'^IAI 



— — mKi{m\\\) 
A 



Finally, substituting this expression into Eq. (C.2), we arrive at Eq. (^) of the main text. 
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